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ABSTRACT 

We study crustal oscillations in magnetars including corrections for a finite Alfven velocity. Our 
crust model uses a new nuclear mass formula that predicts nuclear masses with an accuracy very 
close to that of the Finite Range Droplet Model. This mass model for equilibrium nuclei also includes 
recent developments in the nuclear physics, in particular, shell corrections and an updated neutron- 
drip line. We perturb our crust model to predict axial crust modes and assign them to observed 
giant flare quasi-periodic oscillation (QPO) frequencies from SGR 1806-20. The QPOs associated 
with the fundamental and first harmonic can be used to constrain magnetar masses and radii. We use 
these modes and the phenomcnological equations of state from Steiner et al. to find a magnetar crust 
which reproduces observations of SGR 1806-20. We find magnetar crusts that match observations for 
various magnetic field strengths, entrainment of the free neutron gas in the inner crust, and crust-core 
transition densities. Matching observations with a field-free model we obtain the approximate values 
of M = 1.35 M0 and R = 11.9 km. Matching observations using a model with the surface dipole 
field of SGR 1806-20 {B = 2.4 x lO^^ G) we obtain the approximate values of M = 1.25 Mq and 
R = 12.4 km. Without significant entrainment of the free neutron gas the magnetar requires a larger 
mass and radius to reproduce observations. If the crust-core transition occurs at a lower density the 
magnetar requires a lower mass and a larger radius to reproduce observations. 
Subject headings: dense matter — stars: neutron — stars: magnetic fields — stars: oscillations 



1. INTRODUCTION 

Highly magnetized and isolated neutron stars, known 
as magnetars, emit irregular and extremely energetic 
gamma ray flares. These flares are thought to occur 
following a starquake, in which the magnetar's crust 
is fractured by a reconfiguring magnetic field. Quasi- 
periodic oscillations (Q POs) are observed in the tails 
of giant flare emissions (|Barat et al. 1983 Israel et al. 



2005 



2006 



IStrohmayer &: WaUs||:jO(iSl |Waiis ^ ^irobmayer 
Torsional modes of the oscillating crust suc- 



cessfully explain most QPOs ( Duncan 1998 



Piro 2005; 



Strohmayer fc Watts 2006| Samueisson &: Andersson| 



2007p . Another model associates QPOs with magneto 



hydrqdynamic (MHD) modes in the core ( Glampedakis 
et al.|2006 Levin|2006| . Neither model is able to predict 
ah ot the observed moc 



e frequencies. Crustal oscillations 
cannot easily reproduce all of the low-frequency modes, 
and core MHD modes are una ble to reproduce the high - 
est QPO frequencies observed (van Hoven & Levin 2012 ). 

Resolving these two models and obtaining a complete 
understanding of the mode structure in a highly mag- 
netized neutron star remains challenging. Nevertheless, 
if some of the observed QPOs are associated with the 
torsional modes of the neutron star crust, then magne- 
tars can give a unique insight into the microphysics of 
the neutr on star cr ust, e.g., the nuclear symmetry en- 
ergy (Steiner & Watts 2009). In this work, we assume 
that one ot the lower QPO frequencies and the 626 Hz 
mode observed from SGR 1806-20 both represent the tor- 
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sional modes of the crust. If this is indeed the case 
constraints o n the mass and radius of the magnetar can 
be obtained ( Samuelsson fc Andersson|2007[ [Lattimer fc 



P rakashpOO? 

'I'his study uses modern equations of state (EOSs) to 
predict torsional mode frequencies and compares them to 
observed QPO frequencies to determine magnetar masses 
and radii. The crust EOS is based on a liquid droplet 
model which predicts nuclear masses to within 1.2 MeV 
close to the accuracy obtained in the 



( Steiner 



2012 



i-'mite Range Droplet Model (MoUer et al. 1995). The 



core EOS is based on recent neutron star mass and ra- 
dius constraints from observations of photosphere radius 
expansion bursts (PREs) and the q uiescent emission of 



low-mass X-ray bi naries (LMXBs) (Steiner et al. 2010 



Steiner et al.|2012 ) . At each baryon density we compute 



the shear modulus based on the equilibrium configura- 
tion of nuclei. This study also examines the effects of a 
magnetic field in an oscillating neutron star crust model. 
By a dding the effect of themagnetic field on electrons. 



Broderick et al. 



([2000|, we revise the m agnetic 
composition of the crust! Lai fc Shapiro 19911 with a 
new determination of equilibrium nuclei, 'lb describe 
the oscillating crust we use the axial pert urbation equa- 
tions with general rclati vistic corrections ( Schumaker fc 



Thorne 1983;, Samuelsson fc Andersson 2007) and include 
the magnetic field as done tor the non-relativistic case 



( [Pirol[2005l [steiner fc Watts|[2009l ). 

In §2 we present the magnetized crust composition 
based on our mass model (a detailed description of this 
formalism is in Appendix A). Section 3 contains a sum- 
mary of the axial perturbation equations for the crust 
modes. We then use, in §4, the predicted fundamen- 
tal and harmonic frequencies, along with the magnetized 
crust composition, to constrain the masses and radii of 
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Fig. 1. — Equilibrium composition of the crust without a mag- 
netic field for a model that neglects shell effects (left panel) and 
for a model that includes them (right panel). The blue lower curve 
corresponds to the proton number (Z) and the black upper curve 
corresponds to the neutron number (A'^). 

magnetars. 

2. CRUST COMPOSITION 

For an isolated neutron star we determine the crust 
composition by finding the equilibrium nucleus at a given 
baryon density. At the temperatures of interest, typi- 
cally T < 10^ K, the crust consists of only one nucleus at 
a given depth. The outer crust consists of a lattice of nu- 
clei embedd ed in a degenerate electron Fermi gas (Baym 
et al.|1971a ). The neutron-drip point, the point at which 
it becomes energetically favorable for neutrons to drip 
out of nuclei, defines the boundary between the outer 
and inner crust. The inner crust can then be described 
as a lattice of nuclei embedded in both an electron and 

As described in Ap- 



:as (Baym et al. 1971a 



neutron 

pendix \M at a given proton number Z, atomic number 
A, and baryon density n, the total energy density of the 
crust will have contributions from the nuclear binding 
energy, the Coulomb lattice, the electron gas, and the 
neutron gas. At a given baryon density the equilibrium 
nucleus minimizes the total energy density of the system. 
The most energetically favorable nuclei tend to contain 
a closed shell of protons or neutrons due to shell correc- 
tions. As density increases equilibrium nuclei will move 
to higher closed shells of protons and neutrons with the 
most neutron-rich nuclei seen between the neutron-drip 
point and the crust-core transition. The above features 
of the crust composition can be seen in Figure [l] We 
ignore the deformation of nuclei at high densities in the 
crust composition. 

In the outer crust a strong magnetic field will force 
electrons to occupy the lowest Landau levels. At higher 
baryon densities electrons can occupy higher Landau lev- 
els and thus their energy density approaches the field- 
free case. For this reason, as seen in Table [l] only the 
outer crust equilibrium composition is significantly al- 
tered. For B < 10^* G we can ignore both the effect of 
the magnetic fie ld on the structure of the nuclei in the 



crust (see, e.g., Harding & Lai 2006 Nag et al 
and on the gross structure ot the neutron star 



20091 



ardall 



et al.|[200T ). Qualitatively s imilar results have also been 
obtained in Chamel et al. (2012) using Hartree-Fock- 



Bogoliubov instead of our liquid droplet model. 

3. TORSIONAL OSCILLATIONS IN A STRONG MAGNETIC 

FIELD 

We describe the axial crust mod es of an o scillating 
neutro n sta r following the work of ,Schumaker fc Thorne] 
(19831 and Samuelsson & Andersson (2007r We com- 
bme two forms ot the axial perturbation equation. In 
the non-magnetic case, the equation for the axial per 



turbation $ can be written in the form ([Samuelsson fc 

0, in which primes in- 



Andersson|2007l ) C' + F'^ + 



dicate derivatives with respect to the radial coordinate, 
and F and G are functions of the shear velocity Vs and 
the metric functions ly and A for a static and spherically 



symmetric spacetime metric, ds = 



(d0^ 



sm 



The radial component of the per- 



turbation is suppressed due to the strong vertical strat 
ification. Working in the isotropic limit, we incorporate 
corrections for a finite Alfven velocity va — B/JA-Kpj 



by analogy with the N ewtonian expressions (jPiro 2005 
Steiner fc Watts||2009[ ). The resuh is 



{vl + vDi" + vl^ {In [r^e-^ [e + p) vl] ] C 



+ e 



2A 



2 

e uj 



d 

' dr 



e = o.(i) 



In this expression r is the radius, e is the energy density, 
p is the pressure, uj is the angular frequency, and I is 
the angular wavenumber. The shear velocity is Vg — 
^ [ij p. Here /i is the shear modulus, for which we use 
the for mulation appropriate fo r a body-centered cubic 
lattice ( Strohmayer et al.]|1991 ). 



0.1194r 



i + o.595(ro/r) 



(2) 



We integrate equation ([T]) over the solid crust, which lies 
between the crust-core interface at r = i?core and where 
the lattice melts at r = i?crust- The melting transition is 
determined by where the pla sma coupling parameter F = 
{Zef lak^T = V^^xt = 175 ( [Farouki fc Hamaguchi|l993 



Horowitz et al. 



2010|. Here a = (3/47rni)^^is the radius 



ot the Wigner-seitz cell, Z is the atomic charge number. 
Hi is the ion number density, and the temperature is T 
= 3.0 X 10^ K. 

For the boundary conditions needed to solve equa- 
tion we require the traction, to vanish at the 
top and bottom of the crust. This is a good approxima- 
tion near the surface where pressure vanishes. The de- 
scription of matter near the crust-core transition is com- 
plicated by the appearance of nuclear pasta. Since the 
quasi-free neutrons are superfluid, assuming the traction 
vanishes at the crust-core boundary may also be a good 
approximation. An additional impact of the superfluid 
is that some fraction, /ent, of the quasi-free neu t rons are 
entra ined with the nuclei (Chamel et al. 2012 Chamel 
20131. We assume zero traction at the crust-core transi- 



tion and leave a more complete description of matter at 
the highest densities to future work. 

For a given Z, equation ([T]), when integrated over the 
crust with the boundary conditions described here, has 
an eigenvalue uj that is uniquely determined by the crust 
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TABLE 1 

Magnetic Equilibrium Nuclei Below Crust-Core Transition 



Pmax 



Nuclei 

50rp- 

84o 
340e 

iGe 
iNi 

36 '^^ 

40 



B, = 
3.96 X 10* 
1.10 X 10^ 
1.41 X 10^ 
5.52 X 10^ 
1.69 X lOio 
3.88 X 10^° 
1.52 X 10" 
1.85 X 10" 

2.48 X 10" 

3.49 X 10" 
5.15 X 10" 
1.18 X 10^2 
2.34 X 10^2 
3.81 X 10^2 
8.32 X 10^2 
1.06 X 10" 



B-t = 10 

4.15 X 10* 

1.16 X 10^ 



B* = 10^ 
1.80 X 10^ 
2.79 X 10^ 

4.76 X 10^ 
1.61 X lOi'' 



B* = 10^ 
9.90 X 10^ 
1.96 X 10^0 

5.20 X 10^0 
7.32 X lO^o 
9.81 X 10" 
1.60 X 10" 
1.68 X 10" 
1.94 X 10" 
3.85 X 10" 
5.68 X 10" 



Note. — We adopt the format of Lai & Shapiro (1991) where pmax is the maximum mass density where the equilibrium nucleus is 
present. If a density value is unchanged the following column is blank. Here Bt = 5/(4.414 X 10^'' G), which is the ratio of the magnetic 
field to the critical field, defined as the field at which the cyclotron energy equals the electron rest-mass. 



thickness A — Rr 



■ Rrnrr and thc ncutron star radius. 



These in turn depend on the equation of state. 

4. MAGNETAR MASSES AND RADII 
For the core, we use t he probability distri bution for 



Steiner et al. (2012) from ob- 
rom the quiescent emission of 



the EOS determined by 
servations of PREs and 
LMXBs. In this work, we use a range of five high-density 
EOSs corresponding to the most probable equation of 
state along with its 1 and 2-a lower and upper bounds. 

QPO frequencies have been detected in two magne- 
tars, SGR 1806-20 and SGR 1900-H4 dlsrael et al 



2005] IStrohmayer fc W atts"2005| [Watts fc Strohmayer 
20061 jStrohmayer & Watts 2006"f. The 29 Hz mode in 
SGR 1806-20 and thc 28 Hz mode in SGR 1900-t-14 are 
often assumed to be the fundamental torsional modes, 
but an 18 Hz mode was also observed in SGR 1806—20 
and a lower frequency mode is not ruled out by the 
1900+14 data. SGR 1806-20 also showed a very clear 
626 Hz mode, possibly matching the first radial harmonic 
(n = 1). Several other modes are observed between 50 
and 200 Hz, and these can be matched with higher an- 
gular momentum harmonics, I > 0. 

Each equation of state gives mass-radius combinations 
with different crust thicknesses, and hence a unique fun- 
damental mode {n = 0) and harmonic mode {n = 1). 
We can constrain the masses and radii of magnetars by 
matching predicted fundamental modes and harmonic 
modes to observed QPOs. While GR tends to de- 
crease the frequencies, softer core equations of state with 
smaller radii tend to increase the frequencies. Because of 
this latter effect, we get frequencies which are larger than 



that obtained in Steiner fc Watts] (2009). The n = 
I = 2 mode corresponds to the 29HzlJPD of SGR 1806- 
20. Our model predicts n = 1 harmonic modes near 
600 Hz and we compare these predicted modes with the 
626 Hz QPO of SGR 1806-20. 
To find crusts with fundamentals that match the 29 Hz 



QPO we model crust perturbations in magnetars be- 
tween 0.8-2.0 M0 with magnetic fields matching the sur- 
face dipole field of SGR 1806-20. Whichever crust has an 
n = 0, ^ = 2 mode that matches the 29 Hz QPO we take 
as the crust of the magnetar. This method is demon- 
strated in Figure M where crusts are constructed using 
the SLy4 crust EO^ (jChabanat et al.|1995|). Crusts with 
harmonics that matcnthc 626 Hz QPO are found using 
an identical technique. Whichever crust has a n = 1 
mode which reproduces the observed QPO we take as 
crust of the magnetar. The same analysis is repeated for 
the 1 and 2-a lower and upper bounds on the core EOS. 

A comparison of masses and radii from fundamental 
and harmonic modes can be seen in Figure |3] The inter- 
section of fundamental and harmonic masses and radii 
on the mass versus radius plot gives a crust that best 
matches the properties of SGR 1806-20. The mass and 
radius found for SGR 1806-20 depend on the proper- 
ties of the interior of the magnetar which determine the 
fundamental and harmonic modes. This study focuses 
on varying three aspects of the interior physics that re- 
main unknown, namely, the magnetic field strength in 
the crust, the crust-core transition density, and the de- 
gree of free neutron entrainment in the inner crust. 

First, we examine the sensitivity of fundamental and 
harmonic modes to the strength of the magnetic field. 
Strong magnetic fields melt the outermost boundary 
of the crust (i.e., push the melting point of the one- 
component plasma to higher pressures). Since i?core 
remains fixed and i?crust decreases, a strong magnetic 
field thins the crust (i.e., decreases A). Although a 
strong magnetic field can decrease the crust thickness 
and change the composition of the outer crust, the over- 
all impact on predicted fundamental and harmonic mode 
frequencies is negligible. We find that predicted funda- 
mental modes from magnetized crusts are nearly iden- 
tical to the field-free case. The magnetic field is not 
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Fig. 2. — Frequency of the fundamental I = 2 mode as a function 
of the magnetar mass fo r th e core EOS prob ability distribution 
(centroid and ±2(7) from Stcincr ct al. (2012) and an SLy4 crust 
EOS. The dashed black line indicates the observed 29 Hz QPO 
of SGR 1806-20. The frequencies are evaluated for a crust-core 
transition density of 0.12 fm~^ with U = G. 

a determining factor because fundamental crust modes 
are entirely set by our choice of radius for the magne- 
tar. That is, fundamental modes are entirely set by the 
equation of state. However, a magnetized crust can sig- 
nificantly alter predicted harmonic modes. The n = 1 
modes are sensitive to the magnetic field, especially in th e 



outer crust where va > Vg ( Piro|2005 



Nandi et al. 2012) 



i 



The magnetized crusts that match observed yJ^Us can 
be seen in Figure [3] 

We test two extreme values for the crust-core tran- 
sition density, nt, from Oyamatsu & lida (|2007|). The 



exact value of the crust-core transition density is un- 
known, in part because the density dependence of the 
nuclear symmetry energy in the inner crust is not well 
constrained and also b ecause of the possib le existence of 
nuclear pasta (see e.g. Newton et aL||2013 |. In Figure [s] 
for each transition density we examine various magnetic 
field strengths to find masses and radii from intersections 
of fundamentals and harmonics. For a larger crust-core 
transition at 0.12 fm""^ and at the surface dipole field of 
SGR 1806-2(0 {B = 2.4 x 10^^ G) we find the magnetar 
to have M = 1.25 and R = 12.4 km. We must ex- 
trapolate outside the equation of state curves to approx- 
imate a mass and radius for the lower crust-core transi- 
tion density 0.08 fm""^. This crust-core transition gives 
M = 0.96 M© and R = 13.5 km for SGR 1806-20. In 
either case, if we assume that the magnetic field inside 
the crust is larger than the observed surface field, then a 
smaller mass and larger radius is implied. If the magnetic 
field is too large, the implied radius will be far outside 
radii implied by mass and radius ob servations from th e 
quiescent LMXBs in Ml 3 and a; Gen (iSteiner et al.|2010|). 



4 McGill SGR/AXP Online Catalog, 

[htt p: //www .physics .mcgill . c a/ -p ul sar /magnet ar /main ■ htrnl] 



As can be seen in Figure [3] only harmonic modes are af- 
fected by a change in the crust-core transition density. 
A change in the crust-core transition density will change 
the crust thickness and harmonic modes scale with the 
crust thickness. Fundamental modes remain unchanged, 
however, because they scale with the radius of the entire 
star. 

Entrainment of the free neutron gas in the inner crust 
alters both fundamental and harmonic modes. We deter- 
mine the degree of entrainment by the fraction of the free 
neutron gas that moves with the nuclei lattice during a 
crust oscillation. The intermediate cases of /o„t — 0.75 
and 0.50 can be seen in Figure|4] If a lower fraction of the 
free neutrons are entrained, then larger masses and radii 
are implied for the magnetar. It can be seen in Figure |4] 
that fundamental modes are significantly altered by the 
degree of entrainment. This occurs because the funda- 
mental mode energy is concentrated deeper in the crust, 
whereas harmonic modes have their ener gy distributed 
more uniformly over the crust |Piro||2005| ) . We find that 
a large degree of entrainment, Jent > 0.75, is required for 
predicted crust modes to match observed QPOs using 
crust models built upon the SLy4 interactio n (Chabanat 



et al. 19951). In contrast, for the Rs model (l^'riedrich &; 



Reinhard 



Te 



"03 



986 1 a low degree of entrainment is required 



irnj.3o 



For comparison, we repea t the above analysis with th e 
Rs crust equation of state (Friedrich & Reinhard] 1986 1 . 
This EOS gives frequencies between 15-20 Hz for the 
n — 0, 1 — 2 fundamental mode. Therefore, we must as- 
sign the predicted fundamental modes from Rs to the ob- 
served 18 Hz QPO. The density derivative of the nuclear 
symmetry energy, L, is different in the Rs model and the 
SLy4 model. The Rs model corresponds to L = 86 MeV 
whereas the SLy4 model has L — 46 MeV. The Rs model 
has a smaller fundamental frequency than SLy4 because 
it has a larger L value. 

The Rs EOS also gives harmonic modes consistent with 
the 626 Hz QPO. Using the Rs EOS and a larger crust- 
core transition at 0.12 fm^^ and at the surface dipole field 
of SGR 1806-20 {B ^ 2.4 x 10^^ G) we find the magnetar 
to have M = 1.10 M© and R = 13.8 km. In Figure [5] we 
show the fundamental and harmonic mode intersections 
for various magnetic fields. Finally, we examine the effect 
of entrainment on frequencies predicted by the Rs EOS. 
The predicted n = 0, Z = 2 fundamental gives frequencies 
near the 29 Hz QPO of SGR 1806-20 when there is a low 
degree of entrainment of the free neutrons. The behavior 
of the fundamental with the degree of entrainment can 
be seen in Figure |6] For example, for fent = 0.25 we find 
the magnetar to have M = 1.12 M© and R = 12.0km. 

5. DISCUSSION 

Fundamental and harmonic modes are altered by the 
presence of a magnetic field, by the degree of entrainment 
of the free neutron gas in the inner crust, and by the loca- 
tion of the crust-core transition. Although the observed 
surface dipole field strengths are too weak to affect the 
fundamental torsional modes of magnetars, the harmonic 
modes are significantly altered by fields > 10^^ G. For a 
higher transition density at 0.12 fm~^ a magnetic field of 
B > 4.0 X 10^^ G gives no mass and radius solutions con- 
sistent with the mass and radius constraints from PREs 
and LMXBs. For larger magnetic fields extrapolation 
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Fig. 3. — Magnetar mass as a function of radius for tiie core EOS probability distribution from [Steiner et al.| l |20l"2] |. Frequencies are 
evaluated using the SLy4 crust EOS for two different crust-core transition densities, rat = 0.12 fm~'^ (left plot) and rit = 0.08 fm~^ 
(right plot). The thick red solid line indicates masses and radii for which the fundamental mode has a frequency of 29 Hz. The black 
short-dashed line indicates masses and radii for a 626 Hz harmonic mode and _B = G. Masses and radii from 626 Hz harmonic modes 
with magnetized crusts are labeled accordingly. Arrows indicate masses and radii that match both the fundamental and the harmonic 
modes for the field-free case and the case with the magnetic field of SGR 1806-20 {B = 2.4 X 10^^ G). 
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Fig. 4. — Magnetar mas s as a function o f radius for EOS proba- 
bility distribution from Stci ner et al.||2bl2| . Frequencies are evalu- 
ated using the SLy4 crust EOS with B = G and m = 0.12 fm"^. 
The red dot-dashed, blue dotted, and black dashed lines indicate 
masses and radii from fundamental modes of frequency 29 Hz for 
different free neutron entrainment fractions /ont • The shaded band 
indicates masses and radii from 626 Hz harmonic modes as /end 
is varied from 0.50 to 1.0. Arrows indicate the masses and radii 
that match both the fundamental and the harmonic modes for 
/ont = 1.0, 0.75, and 0.50. 
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Fig. 5. — The same as Fig.js] but for the Rs crust EOS with rat = 
0.12 fm~^ and for a 18 Hz fundamental frequency. The thick red 
solid line indicates masses and radii determined from fundamental 
modes. Masses and radii from harmonic modes with magnetized 
crusts are labeled accordingly. Arrows indicate masses and radii 
that match both the fundamental and the harmonic modes for the 
field-free case and the case with the magnetic field of SGR 1806-20 
{B = 2.4 X 10^5 G). 

outside of the EOS curves is necessary or another EOS is 
needed. The same is true for a lower crust-core transition 
at 0.08 fm"^ and a magnetic field of B > 2.0 x G. 
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Fig. 6.— The same as Fig.|4] but for the Rs crust EOS with B = 
OG. Here the free neutron entrainment fraction /ent is varied from 
0.20 to 0.30, with /ent labeled next to the corresponding curves. 
The red dot-dashed, blue dotted, and black dashed lines indicate 
masses and radii from the 29 Hz fundamental mode. The shaded 
band indicates masses and radii from the 626 Hz harmonic mode. 
Arrows indicate the masses and radii for f^nt = 0.30, 0.25, and 
0.20 that match both the fundamental and harmonic modes. 

For both transition densities we find crusts for tlie field- 
free case and for various magnetic fields as seen in Figure 
[Sj For all transition densities a magnetized crust requires 
a lower mass than the field-free case in order to contain 
a mode consistent with the observed 626 Hz QPO. The 
field-free case gives a minimum radius for a crust that 
can reproduce observations of SGR 1806-20; our model 
requires i? > 11.9 km for SGR 1806-20. 

The sensitivity of harmonic modes to the magnetic 
field strength hints at the behavior of the magnetic field 
with depth. If the value of the magnetic field is constant 
throughout the crust the surface dipole field must be con- 
sidered an upper-limit on the field strength in order to 
find mass and radius solutions consistent with constraints 
from PREs and LMXBs. The surface dipole field must 
be a large overestimate of the surface field if the magnetic 
field were to be larger deeper in the crust. However, the 
model could allow for a larger surface dipole field than is 
observed if the magnetic field were weaker deeper in the 
crust. Alternatively, any configuration of the magnetic 
field is allowed as long as the field strength is < 10^^ G 
everywhere in the crust. 

The degree of entrainment of the free neutron gas in 
the inner crust will alter both fundamental and har- 
monic modes. Crusts with less entrainment require larger 
masses and radii to have modes consistent with observed 
QPOs. We must extrapolate outside the EOS curves to 
approximate a mass and radius for the case of /ont — 0.50 
which requires a mass of M — 1.83 M0 and a radius of 
R = 13.9 km; whereas /ont = 1-0 gives AI — 1.35 Mq 



and R = 11.9 km. The comparison of entrainm ent frac- 
tions can b e seen in Figure [4] It was found in Chamel 
et al. (2012) that entrainment makes global shear modes 
inconsistent with observed QPOs. In contrast, we find 
that at least approximately /ent = 0-75 is required to 
have modes consistent with QPOs. 

The value of the crust-core transition density can sig- 
nificantly alter harmonic modes. A larger transition den- 
sity will increase the crust thickness by decreasing i?coro- 
The thicker crust will have a larger shear velocity on 
average when integrating a perturbation over the crust. 
Fundamental modes are only slightly altered by a high 
transition density. These features can be seen when com- 
paring both panels of Figure [3] 

Magnetar giant flare QPOs give us a unique oppor- 
tunity to probe the neutron star mass-radius relation. 
Fundamental torsional modes are largely independent of 
the crust-core transition density and the magnetic field 
strength. Harmonic modes are sensitive to the details 
of the neutron star interior: the transition density, en- 
trainment of the free neutrons, and the magnetic field 
strength. Comparison of fundamental and harmonic 
modes gives solutions for magnetar masses and radii. 
Since fundamental and harmonic modes are set by differ- 
ent properties of the neutron star, comparison of these 
modes gives us insight into the magnetar interior. In par- 
ticular, we can find the required interior conditions that 
give solutions for magnetar masses and radii that are 
consistent with mass and radius constraints from PREs 
and LMXBs. Solutions most consistent with these con- 
straints have a large degree of free neutron entrainment, 
large crust-core transition densities, and lower magnetic 
fields than observed surface dipole fields. 

Some of the basic trends observed in t he results above 
have also been observe d bylSotani et al. j j2013) . We also 
find, in agreement with Sotani et al. ( |20l3 ), smaller fun- 
damental mode frequencies for crust EUSs with larger 
values of L, as SLy4 corresponds to L = 46 MeV, and Rs 
corresponds to L = 86 MeV. Both works find that more 
entrainmen t decreases the fund amental frequency (In the 
notation of iSotani et all ( |2013l) , Ns/Na = 1 - /ont)- Our 
work includes nuclear snell enects in a more consistent 
fashion, an d thus it is m o re dif ficult to vary L continu- 
ously as in Sotani et al. (2013). Also, we only employ 
EOSs that are consistent with recent constrain ts from 
neutron star mass and radius measurements (Steiner 



et al. 2012) that rule out larger values of L. One criti- 
cal tfimg missing in both works is a complete evaluation 
of how the entrainment in the crust might be correlated 
with L. Work in this direction is in progress. 
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DOE grant bE-FG02-00ER41132, Chandra grant TMl- 
12003X, NSF grant 1109176, and by the Joint Institute 
for Nuclear Astrophysics at MSU under NSF PHY grant 
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APPENDIX 

THE CRUST EQUATION OF STATE 

To compute the energy density of matter in the crust w, we start with an expression similar to that used by |Bayni| 
et al. ( 1971a[ ). We take our crust to be composed of "drops" of nuclear matter with volume fraction Xi within the 
nucieias the density of neutrons and protons are n„ and Hp, respectively, and we denote ni = n„ + np to be the average 
baryon density inside a nucleus. The dripped neutrons, with density ridrip, occupy a fraction 1 — x of the volume. The 
density of nucleons per unit volume is thus n = xiP'n + np) + (1 — x)ndrip, and the density of electrons is rig. As the 
density approaches nuclear saturation the fraction of space filled by the neutron gas approaches unity. 

The energy density w has contributions from nuclei (including the Coulomb lattice contribution), dripped neutrons, 
and electrons: 



w{Z,A,n) = X 



+ npTUp + ni 



EhindjZ, A) 

A 



+ (1 - x)e(n„ 



, Up = 0) + We(ne). 



(Al) 



This expression is valid for any baryon density below the transition density (« 10^'*gcm~^). Here e(n„,np) is the 
energy density, including rest mass, of homogeneous bulk matter at a given neutro n and proton number d ensity. We 
compute e using the bulk matter Hamiltonian in the Skyrme model ( Skyrme||1959 ) with SLy4 coefficients (Chabanat 
eFaL][T995l ). 
I'he energy density of the nucleus is 



«nW„ + UplTlp + ni - 



ad(^, A) 

A 



ein^,np) + '^iE^. 



hcU 



+ £^pair) + WCoul. 



(A2) 



In this expression, n„ and Hp are the neutron and proton densities inside th e nucleus. For the nuclear and lattice 
contributions to the energy density -Ebind, we use a liquid-drop mass model (jBaym et al. 1971a|b Ravenhall et al. 
1983 Steiner| [2008 ) that includes the lattice contribution in the Coulomb term wcoui., as well as surface (i/'surf), shell 
(-fc'sheii)i and pairing (-Epair) corrections to the homogeneous bulk matter Hamiltonian e. At lower densities, the energy 
per particle in the crust is minimized when ridrip = 0, and after the neutron-drip point (about 4 x 10^^ g/cm'^), the 
energy per particle is minimized only when ndrip > 0. The baryon number density inside a nucleus ni is determined 
from 



ni = no + nil 



(A3) 



where / = 1 — 2Z /A is the isospin asymmetry, no is the nuclear saturation density of bulk homogeneous matter, and 
ni < is a correction due to both the isospin asy mmetry, which decreases the saturation density, and the Coulomb 
interaction, which increases the saturation density (Steiner 2008). The average neutron and proton densities within 
the nucleus are then determined from ni and / via 



ni , 
— 1 
2 ^ 



2 ^ 



(A4) 



where rj = 5/1 = 0.92 is a constant of our model that determines the thickness of a neutron skin ( Steiner||2008[ ) , i.e. 
the difference between neutron and proton radii, and 5 = 1 — 2np/{nn + np) is the density asymmetry. 



The next three terms in equation (A2) are the surface, shell, and pairing corrections. The surface correction is 



proportional to the surface tension cr, the nuclear surface area A^^^ , and density asymmetry (5, 



Esmi — o 



(m-nA^ 



1/3 



fl - 0,5" 



(A5) 



where > is a parameter that represents the surfa ce asymmetry ([M yers fc Swi atecki|1969 Steiner et al.|2005 1 . The 
shell correction to the binding energy per baryon is ( jPieperink fc van lsacker^^20DU ) 

(A6) 



where 



S3 = 



Tlyfly Zy Zy 



and 



Zy 



\D, - Zy. 



(A7) 
(A8) 
(A9) 



(AlO) 
(All) 
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TABLE 2 
Parameters of the mass model. 

parameter value 



no 


0.1740 fm"^ 


ni 


-0.0157 fm-3 


V 


0.9208 


OS 


1.964 


a 


1.164 MeV 


ai 


-1.217 MeV 


02 


0.0256 MeV 


£13 


0.0038 MeV 




0.0357 MeV 


dp 


5.277 MeV 



The parameters Dn and Dz correspond to the degeneracy of the neutron and proton shehs, i.e., the difference between 
the magic numbers enclosing the current amount of neutrons or protons. The quantities n„ and Zy are the number of 
valence neutrons and protons, i.e., the difference between the current number of protons or neutrons and the preceding 



magic number. The pairing contribution to the nuclear binding energy is taken from Brehm ( 1989 1 with updated 
coefficients, 

-apA"^/"^, even-even 
impair = { +apA~^/^, odd-odd , (A12) 
0, even-odd 



where Op is a constant of our model. The last term in equation ( A2 ) is the Coulomb energy density. 



2tt 2 



'Rl(2-3x'^' + x), (A13) 



where e is the Coulomb coupling and Rp is the proton radius (3Z = AirnpRp). The respective % terms in parentheses 
correspond to the Coulomb contribution, the lattice contribution, and a correction that accounts for the filling fraction 
X of the nuclei. Table [2] lists the values of the coefficients used in this mass model. 

The electronic contribution to the energy density is that of an electron gas embedded in a uniform magnetic field. 
The electrons acquire an effective mass m^: in the presence of the magnetic field 

m) = ml^2ml(^+]^ + ]^i^ B^, (A14) 

wher e rrie, x, and y are respectivelv the elect ron mass, principal quantum number, and electron spin along the magnetic 
field ( |Rabi||l928| [Ventura fc Potekhin ' 2001 ). Here = heB/mlc' = B/(4.414 x 10" G) is the ratio of the magnetic 



field to the criticaTficld, defined as the field at which the cyclotron energy equals the electron rest-mass. The electron 
number density and energy density are found by summing over electron states and spins in the limit /ig 3> ruf. 



REFERENCES 



Barat, C, et al. 1983, A&A, 126, 400 

Baym, G., Bethe, H. A., & Pethick, C. J. 1971a, Nucl. Phys. A, 
175, 225 

Baym, G., Pethick, C, & Sutherland, P. 1971b, ApJ, 170, 299 
Brehm, J. J. 1989, Introduction to the Structure of Matter 
(Wiley) 

Broderick, A., Prakash, M., & Lattimer, J. M. 2000, ApJ, 537, 351 
Cardall, G. Y., Prakash, M., & Lattimer, J. M. 2001, ApJ, 554, 
322 

Ghabanat, E., Bonchc, P., Haensel, P., Meyer, J., & Schaeffer, R. 

1995, Phys. Scripta, T56, 231 
Ghamel, N. 2013, Phys. Rev. Lett., 110, 011101 
Ghamel, N., Page, D., & Reddy, S. 2012, arXiv:1210.5169 
Ghamel, N., Pavlov, R., Mihailov, L., Velchev, G., Stoyanov, Z., 

et al. 2012, Phys. Rev. G, 86, 055804 
Dieperink, A. E. L., & van Isacker, P. 2009, Eur. Phys. Jour. A, 

42, 269 

Duncan, R. G. 1998, ApJ, 498, L45 

Farouki, R. T., & Hamaguchi, S. 1993, Phys. Rev. E, 47, 4330 
Friedrich, J., & Reinhard, P.-G. 1986, Phys. Rev. G, 33, 335 
Glampedakis, K., Samuelsson, L., & Andersson, N. 2006, 

MNRAS, 371, L74 
Harding, A. K., & Lai, D. 2006, Rep. Prog. Phys., 69, 2631 
Horowitz, G. J., Schneider, A. S., & Berry, D. K. 2010, Phys. 

Rev. Lett., 104, 231101 
Israel, G. L., et al. 2005, ApJ, 628, L53 



Lai, D., & Shapiro, S. L. 1991, ApJ, 383, 745 

Lattimer, J. M., & Prakash, M. 2007, Phys. Rep., 442, 109 

Levin, Y. 2006, MNRAS, 368, L35 

MoUer, P., Nix, J. R., Myers, W. D., & Swiatecki, W. J. 1995, At. 

Data Nucl. Data Tables, 59, 185 
Myers, W. D., & Swiatecki, W. J. 1969, Ann. Phys., 55, 395 
Nag, N., Ghosh, S., & Chakrabarty, S. 2009, Ann. Phys., 324, 499 
Nandi, R., Ghatterjee, D., & Bandyopadhyay, D. 2012, 

arXiv: 1207.3247 
Newton, W., Gearheart, M., & Li, B.-A. 2013, ApJS, 204, 9 
Oyamatsu, K., & lida, K. 2007, Phys. Rev. C, 75, 015801 
Piro, A. L. 2005, ApJ, 634, L153 
Rabi, I. I. 1928, Zeitschrift fur Physik, 49, 507 
Ravenhall, D. G., Pethick, G. J., & Wilson, J. R. 1983, Phys. 

Rev. Lett., 50, 2066 
Samuelsson, L., & Andersson, N. 2007, MNRAS, 374, 256 
Schumaker, B. L., & Thorne, K. S. 1983, MNRAS, 203, 457 
Skyrme, T. 1959, Nucl. Phys. A, 9, 615 
Sotani, H., Nakazato, K., lida, K., & Oyamatsu, K. 2013, 

MNRAS, 428, L21 
Steiner, A. W. 2008, Phys. Rev. G, 77, 035805 
Steiner, A. W. 2012, Phys. Rev. G, 85, 055804 
Steiner, A. W., Lattimer, J. M., & Brown, E. F. 2010, ApJ, 722, 

33 

Steiner, A. W., Lattimer, J. M., & Brown, E. F. 2012, ApJL, 5 



Magnetar Oscillations 



9 



Stcincr, A. W., Prakash, M., Lattimcr, J. M., & Ellis, P. J. 2005, 

Phys. Rep., 411, 325 
Steiner, A. W., & Watts, A. L. 2009, Phys. Rev. Lett., 103, 

181101 

Strohmayer, T., van Horn, H. M., Ogata, S., lyetomi, H., & 

Ichimaru, S. 1991, ApJ, 375, 679 
Strohmayer, T. E., & Watts, A. L. 2005, ApJ, 632, Llll 
— . 2006, ApJ, 653, 593 



van Hovcn, M., & Levin, Y. 2012, MNRAS, 420, 3035 
Ventura, J., & Potckhin, A. 2001, in The Neutron Star - Black 

Hole Connection, ed. C. Kouveliotou, J. Ventura, & E. van den 

Heuvel, 393 

Watts, A. L., & Strohmayer, T. E. 2006, ApJ, 637, L117 



